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Abstract. 

We find that characteristics of quantum tunneling in the presence of chaos can 
be regarded as a manifestation of the Julia set of the complex dynamical system. 
Several numerical evidences for the standard map together with a rigorous statement 
for the Henon map are presented demonstrating that the complex classical paths which 
contribute to the semiclassical propagator are dense in the Julia set. Chaotic tunneling 
can thus be characterized by the transitivity of the dynamics and high density of the 
trajectories on the Julia set. 



Recent studies on tunneling in multidimensions have revealed that the existence of 
chaos affects the signature of quantum tunneling. The observation of purely quantum 
mechanical calculation in chaotic systems shows that tunneling can become chaotic or 
chaos seems to assist tunneling [|], [|, H, El> H, H H> & IE HO' HH- The idea capturing 
such a novel aspect of tunneling looks very attractive, but a direct connection between 
chaos and tunneling can only be accomplished by interpreting the quantum phenomenon 
by the trajectory of the classical dynamics |], |7|, |], pi], [lljj. When one is particularly 
interested in the tunneling process, the use of complex trajectories is essential since the 
transition due to tunneling occurs where the real-valued classical trajectories cannot 
reach. 

The most well-known technique using the complex space is so-called instanton 
method in which the tunneling penetration is evaluated mainly by a single classical 
path moving on the reversed potential [12 , 13| . On the other hand, in chaotic systems, 
it has been found in the time-domain semiclassical analysis that a bunch of complex 
paths almost equally contribute to the tunneling transition between classically forbidden 
regions (10], |TT[. They typically form a tree-like fractal structure in the complex initial 
value plane, and its outstanding appearance compels us to prepare some concept which 
controls dominating complex paths in the semiclassical sum of contributing candidates 
PH |TT| . All the characteristic structures appearing in the tunneling wavefunction in 
chaotic systems originate from it. 



Julia set describes quantum tunneling in the presence of chaos 



2 



However, the Laputa chain, which was so introduced in [10], ITTJ, has been still 



phenomenological so far, and remains even mysterious if no link to some concept 
compatible with the dynamical system theory is made. One may thus naturally ask 
why such a structure play a special role in the complex trajectory description of 
chaotic tunneling, and what sort of mechanism underlies such conspicuous objects in 
the complex plane. The purpose of this Letter is to provide a clear answer to these 
questions. Our final claim is simple and would be rather natural: Julia set is the origin 
of chaotic tunneling. 

Let us begin with introducing the model system we are concerned with. The system 
we study here is a family of two-dimensional area preserving maps, in which the mixed 
phase space being realized in a certain range of the parameter space. The time evolution 
of the phase point (p, 9) is given as the mapping rule as 

( Pn+1 , e n+1 ) = F( Pn , e n ) = (H'( Pn ) - v'(e n ), e n + H\ Pn ) - v'(e n )). (1) 

Here, H(p) = p 2 /2 and V(6) = Ksin8 are the most standard choice, but suitable 
modification or replacement of the kinetic or the potential term is sometimes helpful 
and will be made according to the target of the analysis. 

Since the map model does not have the energy as the Hamiltonian flow problem 
does so, one cannot consider the tunneling through the energetic barrier, which may 
be a normal setting of the tunneling problem. Instead, dynamical confinement due 
to classically disconnected components such as KAM tori and chaotic components in 
the phase space plays the role of barriers, and the quantum transition between such 
invariant regions is regarded as tunneling |TJJ . 

At least in the first setting of the problem, it is not at all obvious that several 
different situations, such as the tunneling transition out of the quasiperiodic region 
into some chaotic component, or its reverse process, or that between different chaotic 
components, could be treated on the same footing. However, as will be described below 
and also become one of the most important point in the present report, the choice of 
initial and final states does not matter to the whole story. 

Typical quantum mechanical wavefunctions in the mixed phase space are displayed 
in Fig. 1. In both models, the tails of the wavefunctions do not monotonically decay 
even in the tunneling regime, rather there appear several unexpected structures; the 
crossovers of the slope, the plateau regions and irregular interference patterns on it. 



All these characteristics are only qualitatively featured |T(J, jlTf , but they are commonly 
observed not only in the dynamical tunneling problem, but also in the energetic barrier 
tunneling [15|]. 

The semiclassical approach, which is extensively developed in recent studies of 
quantum chaos or quantum chaology |16[, works quite well even when one employs it as 
a tool describing the tunneling process. Apart from an added technical (but sometimes 
crucial) problem originating from the Stokes phenomenon frTfl , which we do not enter 
into details here, our task in the semiclassical analysis is essentially the same as the real 
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one, that is, to evaluate the Van-Vleck propagator: 

^n(P0,Pn) « 4»(PD> #o) eX P { -iS n (po, 9 ) \ , (2) 

p =a in J 

Pn=P 

where the summation is taken for all (po, 9q) which satisfy the boundary conditions for 
the initial momentum po = a and the final momentum p n = {3. Here, S n (po,0o) = 
Y%=i[H(9j) — V(9j) + 9j(pj — Pj+i)} is the action along a classical trajectory, and 
A n (p ,9 ) = [2%h(dp n /d9 ) po }~^ represents the amplitude factor associated with its 
stability. 




p p 



Figure 1. Quantum(bold line) and semiclassical wavefunction (dashed line) for (a) 
the model with H Q (p) = t (wWFFi + and V ^ = Ksin9 > where p d = 5, uj = 2 
and K = 1.2, and (b) the model with H (p) = ^ and V(6) = If sin 0, where K = 1.5. 
In both cases, the initial wavepacket is set as ^(p) = S{p). The real- valued classical 
orbits cannot reach the region outside the dashed lines within the time step taken here 
(n = 6 for (a) and n = 5 for (b)). The semiclassical wavefunction is shifted in order to 
clarify the structure. 



Since we here take the ^-representation, p should be a real quantity. So, the 
canonical partner #0 may be used to identify the (complexified) trajectories contributing 
to the sum (0), and it is then allowed to be complex as 9 = £ + ir) (£, r] real). We 
visualize the contributing complex paths by displaying the set ||T0| , [11] : 



M n = |J M*/ = |J {(p,9) e C 2 I Vn = (3} (3) 
pen. /3eR 

on the ^Q-plane of the slice {po = a} for some initial condition a G R. The set Ai„ on 
the #o-plane, which usually looks like clouds or wisteria trellis on a macroscopic scale, 



is decomposed into finer and finer structures as it is magnified [|10| , [11]]. One can see 
that its basic element is a string with various scales. Each string represents a trajectory 
appearing in the semiclassical sum (El). We note that in the integrable limit only the 
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branches connected with the real plane (i.e., rj = 0) survive and all other complicated 
objects disappear |T(], [TT|| . 



A huge amount of candidate paths may discourage us since it appears to be no 
more possible to establish a simple view of tunneling in the presence of chaos. However, 
among all possible candidates the complex paths forming a sequential structure, which 
runs in the vertical direction at the center of Fig. 2(a) and clearly discernible from 
the other aggregated strings, exceed any other candidate paths in amplitude. We have 
called such a characteristic structure the Laputa chain fTO], [□]]. As shown in Fig. 1, one 
finds that semiclassical sum including only such complex paths contained in the Laputa 
chains has reproduced almost all details of tunneling into chaotic regions. Our task is, 
therefore, reduced to clarifying what this marked structure appearing in the initial value 
plane represents. 





Figure 2. (a) A magnification of the initial value representation Ai n on the 0Q-p\a,ne 
in case of the standard map with K = 1.5 and n — 40. The range shown above is 
given as 4.197487346 x 1CT 2 < £ < 4.197487362 x 10" 2 , and 6.693592 x 1CT 4 < 7} < 
6.693601 x 10~ 4 . The initial momentum is set as po = 0. Ai n -set consists of a bunch 
of self-similar objects, whose basic element is a string. Except for the central part, 
such strings are so densely aggregated that individual string cannot be resolved in this 
scale. However, magnifying the black area, one can find a similar structure in the scale 
shown here, that is, the black area is also composed of a bunch of string objects. Each 
string represents an individual component of the semiclassical sum (2). The strings 
running in the vertical direction look as if they cross with each other, but actually 
they avoid with very narrow gaps, that is, the strings form a serial chain-like structure 
connected via narrow gaps, (b) The slice of K + by {po = a}. This is numerically 
obtained by plotting the initial points whose trajectories remain a ball in C 2 with a 
certain finite radius, r = 10 3 , in this case. 



The reason why some complex paths dominate the others in the semiclassical sum 
(0) is, in general, that the imaginary part of their action, lmS n (po,9o), are relatively 
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small. This is because the absolute value of each term in (|2]) is mainly governed by 
Im S n (p , 9 ), rather than the amplitude factor A n (p ,9 ). This in turn means that 
the complex paths forming the Laputa chain should gain small imaginary action as 
compared to the other paths not forming the chain structure. Indeed, as shown later, 
the trajectories initially placed on the Laputa chain approach the real (p, #)-plane 
exponentially, which provide minimal or relatively small imaginary action. 

Conversely, one can say that this property characterizes the Laputa chain and 
makes them distinguishable from the others. Furthermore, they are specific in that 
those trajectories stay in bounded regions because after coming close to the real plane 
they almost follow the behavior of the trajectories on the real plane and the real orbits 
are all bounded in the present situation. 

This is a hint to link the Laputa chain to a proper object compatible with the 
theory of dynamical systems, since the Julia set, which plays a central role in the 
complex dynamical systems, is specified as the set satisfying such a property. More 
precisely, the forward Julia set J + is defined as the boundary of the set K + of points 



whose forward orbits remain in a finite region |18|j : 

K+ = { (p, 9) | {F n (p, 9)} n>0 is bounded } (4) 

and 

J+ = dK + . (5) 

The polynomial diffeomorphism like Henon map /, which is defined on C 2 , has a 
polynomial inverse, so both the forward and the backward iterations can be considered. 
In such a case we define K + (resp. K~) as the set of points in C 2 whose forward (resp. 
backward) orbits are bounded, and J + (resp. J~) to be the boundary of K + (resp. K~) 
which we call the forward (resp. backward) Julia set. The set J = J + n J~ is called the 
Julia set of f. The forward (or backward) Julia set is where the orbits have sensitive 
dependence on initial conditions, which means that the chaotic motion is realized on it. 

Remarkably enough, such a purely mathematical object enters into physics as 
quantum tunneling in chaos. Indeed, as shown in Fig. 2(b), the similarity between 
the chain-shaped structure demonstrated in Fig. 2(a) and the slice of K + by the same 
plane is obvious. The slice of K + shows a typical dendrite-like structure which often 



appears in the one-dimensional complex dynamical systems [19 . The number of strings 
constituting the Laputa chain increases in an exponential rate, so coincidence between 
them becomes better as the time proceeds jHJ ^3|. Note that the orbits put on the 
highly aggregated branches surrounding the Laputa chain do not stay in a finite phase 
space domain but rapidly escape to infinity. 

It is possible to provide a rigorous statement if one focuses on the cubic potential 
model given by putting H(p) = p 2 /2 and V{9) = c9 — 9 3 /3. The map ([I]) is transformed 
to a standard form of the Henon map: 

f:(x,y)\ — ► (y, y 2 + (1 - c) - x), (6) 
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by the affine change of coordinate (p,8) = (y — x, y — 1). The Henon map is known 
to be one of the simplest nonlinear systems in the two-dimensional space, and its 
dynamics is extensively studied by several authors. Among them, investigation from 
the complex dynamical point of view has been developed in the last decade (see, for 
example, |20|, |2l], ^2] and the references therein) by using the pluripotential theory, the 
theory of currents, etc. 

As in the case of the standard map, it is reasonable to focus on the Im^^o,^) 
of each trajectory, but to be compatible with the invariant set of the dynamical system 
one should consider the set of trajectories having the property described above in the 
limit of n going to infinity. The most natural condition would be to select the complex 
orbits whose Im5 n (po, 9 ) has a finite limit even when n goes to infinity. Such a filtering 
only serves as a necessary condition for semiclassically contributing orbits, but it is at 
least true that the trajectories whose Im S n (j>o, 6q) are divergent cannot contribute to 
the semiclassical summation since those orbits either tends to zero in their magnitude or 



will be removed by the Stokes phenomenon [17j . Therefore we define the Laputa chains 
as, 

C L a P uta = {ip, 0) e Moo I Im S n (p, 0) converges absolutely at ip, 9)} (7) 

In this definition M.^ is an object introduced to represent the limit of .M n -set when n 
goes to infinity. More precisely, 

-Moo = (J Ml, (8) 
/3eR 

where M.^ is given as the Hausdorff limit of M*^ = {ip, 6) G C 2 | p n = [3} (compare 
eq. (3)). Thus, the set M.00 corresponds to M. n for the time step l n = 00'. It is possible 



to prove that this Hausdorff limit itself contains the forward Julia set J + [23|, [TlJ, which 
in itself is a partial verification of our numerical observation. The following assertion 
concerning the relation between laputa thus defined and J + is proved by the second- 
named author: 

Theorem. Let F be the time-one map on C 2 associate to the kicked rotor with 
Hip) = p 2 /2 and Vi6) = c9 — 9 3 /3, and h top iF) be the topological entropy of F, 



(i) If h Xap {F\ B? ) > 0, then C Laputa D J + . 

(ii) If F is hyperbolic on J and h top (F\ R 2) > 0, then CLaputa = J + ■ 

(iii) If F is hyperbolic on J and /itop(-F|R 2 ) = log2 ; then laputa = J + ■ 

Here X indicates the closure of the set X. The rough sketch of the proof is as 
follows [^3|, PH : that /ito P (-^|R 2 ) > implies the existence of a saddle periodic point 



X in the real phase space. A principal tool we will employ is the following result which 



is established by Bedford and Smillie f20|, 0, : For a complex one-dimensional locally 



closed submanifold M in either J + or an algebraic variety, there is a constant c > so 
that 

lim — \r n M] = c -dd c G + (9) 

n^+00 2 n 
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in the sense of current, where [M] is the current of integration of M, i.e. [M](0) = 
J m 4>\m, and dd c is the complex Laplacian. In this statement, G + represents the Green 
function for K + given by 



G + {x,y)=\\m^\og + \\r\x,y)\\. 



(10) 



It is easily shown that the support of dd c G + coincides with J + . From this theorem we see 
that the stable manifold of any periodic saddle p is dense in J + , that is, (W s (p)) = J + . 
Using this result, together with the fact that the Hausdorff limit M.^ contains J + , we 
obtain the desired claim. 

This claim gives a mathematical verification to the observation numerically found. 
Indeed, as shown in Fig. 3(a), the trajectories leaving the Laputa chains approach 
exponentially to the real (p, 0)-plane. This makes Im5 n (po,^o) converge absolutely. As 
a demonstration of the theorem, we show in Figs. 3(b) and (c) the set \Ja<8 -M.^ f° r a 
fixed 0o and the slice of J + by {p = a} for the Henon map. One can see how Ai n -set 
shrinks to the J + as a function of n. 
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Figure 3. (a) The distance of from the real plane as a function of time is displayed for 
the orbits whose initial conditions are put on the Laputa chains. The solid line denote 
the case of the standard map. The dotted and broken lines are the ones for the Henon 
map. When the Julia set J exists only on the real plane, the orbits always approach 
the real plane directly (dotted line), but otherwise the orbits first move around in C 2 
space, and then approach the real plane (broken line), (b)-(c) The set U/3</3 •^n /? 
(/?o = 10 10 ) for the Henon map is shown as the solid curves in case of (b) n = 9 and 
(c) n = 10. The slice of J + by {p = a} is shown as the dots in each figure. 



Notice that the assumption h top (F\ R 2) > in the above theorem is a mathematical 
expression which corresponds to the fact that the underlying classical dynamics F\ R 2 
is chaotic. We also note that the slice of the forward Julia set J + by {p = a} can be 
shown to have positive capacity for any initial condition a G R. Thus, the theorem 
above suggests that, unlike the instanton solutions in the integrable bunch of 
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paths in C 2 contributes to the tunneling phenomena if the underlying classical mechanics 
is chaotic. 

It should be noted that the assumption in (i) covers the system with mixed 
phase space, which is the most generic situation in physics. In addition, the physical 
implication or interpretation of another theorem of Bedford and Smillie on the 
transitivity of the dynamics |S0L ETI, [321] is suggestive in our problem. It states that 
for any C 2 -neighborhoods of any two points in the chaotic regions there is an orbit 
in C 2 connecting them, even in the case where the chaotic regions in the real plane 
are mutually disjointed by KAM tori. This property exactly guarantees the transition 
between any disconnected regions on the real-valued classical dynamics, and non-zero 
tunneling amplitude of the wavefunction in arbitrary regions is always realized due to 
the transitivity on the Julia set. 

In this way, with the help of strong mathematical claims, which could be established 
only by extending the dynamics to the complex space, we can clearly understand the 
reason why chaos seems to assist tunneling and can become chaotic; these can be 
attributed to such high density of the tunneling paths in J + and the transitivity of the 
complexified dynamics. So far, the structure of the Julia set has been an object which 
mainly attracts the interest of mathematicians. But the present result implies that 
the Julia set is really observable as chaotic tunneling in various physical and chemical 
phenomena. 
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